Dissipative dynamics of a Harmonic Oscillator : A non-perturbative approach 
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Starting from a microscopic theory, we derive a master equation for a harmonic oscillator coupled 
to a bath of non-interacting oscillators. We follow a non-perturbative approach, proposed earlier 
by us for the free Brownian particle. The diffusion constants are calculated analytically and the 
positivity of the Master Equation is shown to hold above a critical temperature. We compare 
the long time behaviour of the average kinetic and potential energies with known thermodynamic 
results. In the limit of vainishing oscillator frequency of the system, we recover the results of the 
free Brownian particle. 

PACS numbers: To be appear here 



I. INTRODUCTION 

The problem of quantum dissipation has been a long 
standing and important problem and has been stud- 
ied in the context of a variety of phenomena, includ- 
ing relaxation phenomenaP, quantum cosmological 
models fsl, H and in more recent times, in the context of 
cold atomslSj . In contrast to the classical dissipative phe- 
nomenon, the quantum analogue has a number of addi- 
tional constraints, which on the one hand make the prob- 
lem much more complex, and on the other hand makes 
it more interesting. Of these, the first and the foremost 
constraint is that of the non-commutative nature of the 
canonical variables, which leads to stringent conditions 
on the fluctuation of the canonical co-ordinates, which 
is commonly known as the Heisenberg uncertainity rela- 
tions. The second major difference between the classical 
and quantum dissipative dynaimcs, is the fact that the 
autocorrelation fucntion of the fluctuating force has a 
memory effect, which depends on the temperature and 
the Planck Constant h. Thus, while the force autocorre- 
lation becomes delta-correlated (Markovian) at high tem- 
peratures or fi, ^ 0, leading to recovery of the classical 
results, the correlation can have a power law dependence 
on time at very low temperatures 0. Hence, a consis- 
tent theory has to extrapolate between these two extreme 
regimes. Another major issue of the quantum dissipative 
system is the positivity of the density matrix during the 
non-equilibrium dynamics of the system. The positiv- 
ity is necessary for the evolution of the system through 
physical states. 

Historically, several approaches have been taken to de- 
scribe the quantum Brownian Motion 0, [I]- SenitzkyQ, 
applied the quantum dissipation model to radiating elec- 
trmagnetic waves in a cavity and showed that even with- 
out the details of the dissipative medium being known, 
one can formulate a mechanism which leads to quan- 
tum dissipation. One of the first attempts in solving the 
problem was by Zwanzigp^. who considered the long 
range memory effect in the fluctuating force autocorre- 
lation function. Ford, Kac and Mazurpjj] proposed the 
system-bath model for the classical and quantum dissi- 
pative system, where the system is linearly coupled to a 



bath of harmonic oscillators. A formal treatment of the 
system plus bath model, was introduced by Feyman and 
Vernon jl2l]. who derived the influence functional for the 
effective action of the system, by integrating out the bath 
degrees of freedom, using the path integral technique. 
This method was later used by Caldeira and Leggettp^ 
to derive the high temperature behaviour of the master 
equation for a quantum harmonic oscillator coupled to 
a bath of harmonic oscillators. Dekkerfl^ proposed a 
phenomenological master equation where there can be 
diffusion not only in the position space but also in the 
momentum space. Dekker and Valsakumar [Tsj showed 
that the various diffusion constants must obey a cretain 
criterion in order to satisfy the positivity condition as 
shown by Lindbladp^ \v^ . In the weak coupling limit, 
Agarwal[l3 derived some exact results for a dissipative 
harmonic oscillator. Hu, Paz and ZhangflQ^ did an exact 
treatment of the same problem and found that the dif- 
fusion constants are time dependent. In all the previous 
approaches, one of the crucial assumption was that the 
system and bath are initially decoupled and they get cou- 
pled at time t — Q. Hakim and Ambegaokar[20| showed 
that the results would be dramatically different if this 
assumption were not true and that the initial transients 
can play a crucial role in the evolution of the density ma- 
trix. There have been some attempts to derive the master 
equation, phenomenologically, within Lindblad form[2l|. 

In this paper, we derive the master equation of a quan- 
tum harmonic oscillator, coupled to a heat bath, within 
a scheme, proposed earlier by ust22|]. Unlike some of the 
earlier approaches HI, [2^, where the high temperature 
expansion was used to convert the non-local memory ker- 
nel to a local one, we retain the original expression for the 
memory kernel and Taylor expand the non-local terms in 
the local time scheme. Thus, our results can be valid, 
in general, for a wide range of temperatures. The paper 
is organized as follows. In Section II, we formulate the 
problem of the Harmonic oscillator and show how one can 
derive the master equation following our scheme. In Sec- 
tion III, we calculate the diffusion constants and remark 
on their behaviour at high temperatures. We also cal- 
culate the limit in which the positivity condition breaks 
down. In Section IV, we recover the free particle results 
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exactly, in the limit of vanishing oscillator frequency. We 
conclude in Sectilon V. 



II. DERIVATION OF THE MASTER 
EQUATION FOR THE HARMONIC 
OSCILLATOR COUPLED TO A HEAT BATH 

We start with a microscopic model for a harmonic os- 
cillator, which is linearly coupled to a bath of oscilla- 
tors. The total Hamiltonian TC of the system consists 
of three parts, viz., the Hamiltonian for the system Ha, 
the Hamiltonian for the bath Hb, and the interaction 
Hamiltonian 7i/, which describes the linear coupling of 
the system to the heat bath. The three parts of the total 
Hamiltonian Ti. are defined as : 



Ha 
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2M ^ 2 
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where, q and p are the position and momentum coor- 
dinates of the system, with mass M, oscillator frequency 
ujQ, and Qi and Pi are the coordinates of the bath oscil- 
lators, with mass rrii and the frequency coi. The coupling 
constant between the system and the i"* bath oscillators 
is denoted by C;. 

In order to derive a master equation of the reduced 
density matrix, we use the Feynman- Vernon procedure 
to calculate the effective action, after integrating out the 
bath degrees of freedom. We assume that at i = both 
the system and the bath are uncorrelated with each other, 
so that we can write the total initial density matrix of 
the whole system as the direct product of the initial den- 
sity matrices of the system and the bath, respectively. 
Following the work of Feynman and Vernon [l2|, the time 
evolution of the density matrix is given in terms of the 
Influence Functional, as, 

P{<ll,<l2,t) = / dq[dq'2J{qi,q2,t;q[,q2,0)p{q[,q2,0) (2) 



Where the quantity J{qi,q2,t;q[,q'2,0) propagates the 
density matrix p{q'i ,q'2,0) from an initial time to the den- 
sity matrix p{qi, q2, t) at final time, and is given by, 

J{Qi,q2,t;q'i,q2,0) ^ j j 2?gi2?92 exp Qs'eff[gi, 92] 

(3) 

Integrating out the bath degrees of freedom, we obtain 
the nonlocal effective action for the system[l3|, given by, 

^S,s = ^{SA[qi] - SA[q2]) + (Sfl + ^i)dr (4) 



where, 5*^ corresponds to the action of system corre- 
sponding to Eq. [Ta| and, 

1 r 

S_R = --y [q-{T)aR{T - s)q^{s)]ds (5a) 



[q-{T-)ai{T~ s)q+{s)]ds 



(5b) 



here we define, qj- — qi ± q2 and the respective memory 
kernels are given by the relations. 



2rauji 



2kBT 



cos{ujiT) (6a) 



(6b) 



The summations in Eq. [Slcan be replaced by integrals, 
by defining a density of states F{u}) of the bath oscilla- 
tors. The density of states, is defined using the Drude 
form of the cutoffQ, as, 



(7) 



and 



F{uu)C^{Lu)l{2mLu^) = {2M-i /n){u:l/ {lu^ + ojD). (8) 
The memory kernels can now be evaluated as: 
aii{T) = M7cj^ [cot(x) exp(-WcT) 
n-n/x 



1 



exp [- (riTr/x) Wct] 



(9) 



where, x = fvjJc/2kBT . The above expression shows 
clearly that the system has two distinct time scales, one 
given by 1/lOc and the other by h/kBT. We now proceed 
to use our scheme of using a Taylor series expansion, as 
was done in an earlier work|23|. 



00 (;) 



n 



[s-rY 



(10) 



where g^'^ is the Z*'* derivative of the position coordinates 
of the system with respect to time. Inserting Eq. [9] and 
Eq. [To] into Eq. EJ neglecting the total derivative terms, 
we obtain, 



1=0 " 



f''aR{f)df (11) 



For T ^ h/kBT, neglecting the transient terms, we ob- 
tain. 



Sfl(r) 



hx 



•= m=0 



(12) 
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where C is the Riemann Zeta function, and 



and 



S/(r) 



q-{T)q+{T) 



(13) 



where the higher order terms in S/ are negelected as 
they fall off as I/uJc- Note here, that while Eq. [T2l and 
Eq. I13( are local in nature, the nonlocality is implicitly 
present due to the infinite summation and the presence 
of all higher order derivatives of the dynamical variables. 
We now proceed to find a systematic replacement for the 
higher order derivatives. In order to do this we resort 
to the dynamical equations of motion. The dynamical 
equation of motion for the harmonic oscillator can be 
written as. 



where, 



Q 

M 



Q = MQ 



(14) 



l/M \ 



(15) 



The higher order derivatives of the position variable q, 
can be formally obtained from the relation. 



q('+i) = x'Q 



(16) 



Using the eigenvalues and eigenfunctions of the ma- 
trix AA, and performing a similarity transformation, the 
higher order derivatives of the position variable can be 
written as, 



(I) ^ (A2A'x - AiA^) 
^ (A2-A1) 



A2-A1) 



where. 



Ai,2 = -7 ± V7 



- Wn 



(17) 



(18) 



are the eigenvalues of the matrix M. It is interesting to 
note that all the higher order derivatives can be written 
in terms of the canonical coordinates at time t. Inserting 
Eq. [17] into Eq. [121 gives. 



Sfl(r) 



-2kBT-fM 



WAt) + a e_{T)] (19) 



where 



a 



A1A2 
A2 — Ai 



A? 



^coth(^) 



+ 



Ai A2 



^2 



■ coth( 



h2X coth(^ 

1 + 



(20) 



(A2-Al)2 



¥coth(^) 



1 



■ coth(- 



h2X coth(^ 



1 



A? 



(21) 



In the limit of cl;^ ^ oo, a and a' are given by, 

/ X^ An \ 

1 



a 



A2 — Ai 



1 

A? 



^^^coth(,^^^ 



2fcBr 



A1A2 
1 

A| 



ftyAiAj ,ft\/AiA2 
coth(— — — ) 



hX2 



2kBT 



2kBT 



and 



(22) 



A2 — Ai 



2kBT 



coth( 



2kBT' 



WA1A2 ^, ,h^/\l\2. 

COth(— r-; —) — 1 



2fcsT 

?lA2 



2kBT 



2kBT 



coth( 



JlA2_, 

2kBT' 



1 



(23) 



Thus the effective action is given by 
'M . . 



dT 



2kBT-tM 
l9 



-q+q- 



dr 



MujI 



q+q- 



aq^ 



a'ql' 



(24) 



It is worth noting that using our scheme of summing 
the memory effects over the classical paths, we have con- 
verted the non-local action into a local effective action. 
This is the central result of this paper. 

Once the effective action in a local form has been ob- 
tained, viz. Eq. 1241 we follow the prescription given by 
Caldeira and LeggettfTsj. to derive the master equation. 
We calculate the change in the reduced density matrix, 
within a small time interval e, using the functional in- 
tegral method. Collecting the terms which are of the 
order of e, and neglecting higher order terms, we obtain 
the master equation, describing the time evolution of the 
density matrix p{x, y, t) in the following form : 



dp 

dt 



ih 
2M 



- 7(x-y) 



11 



dp dp 

dx dy 
2 



2i 



— + —Xp- ^(x 
dx dy J h? 



- v) 



dp 
dx 



dp 

dy_ 



(25) 
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where the diffusion constants are given by, 



Dr, 



2kTja' 
M 



Dpp = 2kTM-f[a + 47^0;') 



(26a) 
(26b) 
(26c) 



The above form of the master equation is different from 
the one obtained by Caldeira and Leggettfisj. due to the 
fact that our effective action has an extra term depending 
on . It is relevant to note that the form of the mas- 
ter equation, which we have obtained from a microscopic 
model, is similar to that obtained by Dekker 1J|, from a 
phenomenological point of view. However, the functional 
form of the diffusion constants are different from those 
proposed by Dekker. 



III. DISCUSSIONS OF RESULTS 

In this section, we examine the behaviour of the dif- 
fusion constants with respect to temperature as well as 
on the damping and the oscillator frequency. There are 
two dimensionless parameters in this problems, viz., 7/^0 
and kBT/hiUQ. The cutoff frequency lUc is considered to 
be much larger than any of the other scales in the prob- 
lem. In terms of the dimensionless parameter 7/wo there 
are two regimes : (1) the overdamped regime {j/ujo > 1) 
and the underdamped regime (7/(^9 < 1), which shows 
oscillatory behaviour. At high temperatures, the diffu- 
sion constants are the same for both the regimes and are 
given by. 



Dr, 



2kBTMj l^l 

GMkeT 
ZkeT 



^272 



(27a) 
(27b) 
(27c) 



The high temperature behaviour is completely indepen- 
dent of the oscillator frequency, upto the order of l/T. 
In the high temperature limit, the classical diffusive be- 
haviour is recovered, with the well known diffusion con- 
stant, Dpp = iM^ksT , and the anomalous diffusion con- 
stants Dqq and Dpq vanish as l/T. 

The two central issues in the problem of quantum 
brownian motion as discussed earlier are : (1) the Heisen- 
berg Uncertainity Principle and (2) the positivity of the 
density matrix during the time evolution. Dekker and 
Valsakumar[l^ addressed the first issue and have shown 
that the diffusion constants have to satisfy a certain con- 
dition, known as the positivity condition, in order to 
maintain the uncertainity relation. The second issue has 
been addressed by Lindblad[l^ from a more mathemat- 
ical approach, where a general form of the master equa- 
tion was proposed, in order to maintain the positivity of 
the density matrix for all times. Interestingly, if the first 



condition is satisfied, then it automatically guarantees 
the Dekker equation in a Lindblad form[17j|. In terms of 
the diffusion constants the Dekker- Valsakumar positivity 
condition is given bvp^: 



/\ — Dpp Dqq 



Di 



^^72/4 > 



(28) 



At high temperatures, although the anomalous diffusion 
constants vanish as l/T, they are crucial for restoring the 
positivity condition and A approaches a constant value 
^272/12. Using Eq. US] in Eq. ^ we find that the 
positivity condition is only satisfied above a breakdown 
temperature which depends on^ /loq. In Fig. [T]we plot 




FIG. 1: The breakdown temperature Tc as a function of the 
dimensionless parameter 0)0/7- 

the dimensionless quantity ksTc/hj as a function of the 
dimensionless parameter cjq/t- For small wq, this critical 
temperature is of the order of OAh'j /ItB, which agrees 
with the results for the free particle[24|. . From the 




FIG. 2: The average potential energy as a function of the 
temperature for (a) 7/0)0 = 0.01 and (b) 7/00 = 2. The 
solid line is obtained from Eq. [31] and the circles denote the 
thermodynamic result taken from [l[ 

master equation (Eq. I25p , one can also derive the average 
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kinetic and potential energies as the system approaches g^, and qp + pq can be calculated as, 
equilibrium [15(1 . The time derivative of the average values 



j^{Q^) = ^{qP + Pq)+^D,g (29a) 
l(p2) = -MiuKqp + pq) - 47(p2) + 2Dpp (29b) 

j^iqp + pq) = -^(p') - 2Mc^o'('z'> - 27(9^ + pq) - ^d^^ (29c) 



The solutions of the above coupled equation can then be obtained as, 



{q 



7(7 + f])e-^('^-")* ^(^_»)e-2(7+o)« 
^'2M^-^' 2M^^ 2M^^ 



27 27 

- e-2(^-")*7(7 - f^) - C^sf e-2(^+")*7(7 + ^) (30b) 
27 27 

qp + pq) = -2MDgg + Cie-27t + c-ae-^l^-f^)* + ^36-2(7+^^)* (30c) 



Where C^'s are the constants of integration, which de- 
pend on the initial values of the aver age quan tities, viz. 
(g^)o, (p^)o and {qp + pq)o a.nd,n = 1/7^ - uj^. For finite 
oscillator frequency loq, the equilibrium values of the av- 
erage kinetic and potential energies can be obtained from. 



27 

2M27u;2 



■qq) 



(31a) 
(31b) 



1 1 1 1 1 1 1 y| 1 
/ — Dynamic Route 


1 1 1 1 1 1 1 


• Thermodynamic Route 
. 1 . 1 . 1 . 1 . 


1 . 1 . 1 . 1 



At high temperatures, the average value of both ki- 
netic energy {p^)/2M and potential energy Mj^^q^) 12 
approach to ■^ksT, which is in agreement with equipar- 
tition theorem. Thermodynamic quantities of the sys- 
tem can also be calculated from the partition function, 
which can be evaluated by well known technique of imag- 
inary time path integral method^!,, J.]. Above the critical 
temperature (below which the positivity condition is vio- 
lated), equihbrium values of (p^) and (q^) obtained from 
the master equation are in good agreement with those 
obtained from the partition function method as shown in 
Figs. HandO 



FIG. 3: The average kinetic energy as a function of the tem- 
perature for (a) -y/ujo ~ 0.01 and (b) 7/1^^0 = 2. The solid line 
is obtained from Eq. [31] and the circles denote the thermody- 
namic result taken from [ij. 



IV. REVISITING THE FREE PARTICLE : THE 
LIMIT OF VANISHING OSCILLATOR 
FREQUENCY 

It is obvious from the Hamiltonian Eq. [1] that, if we 
take luq ~^ 0, we would end up with the system of a free 
particle connected to a heat bath. This problem was dis- 
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cussed in details by us in an earlier paper. We will present 
here the relevant results for the diffusion constants and 
the averages of and from the earlier paper. The 
diffusion constants are given bv[2^. 



and, 



Daa = 



Dpq — 



M 

2kBTM'y{l + A'y^a'o) 



where, 



coth( 



ksT' 



(32a) 

(32b) 
(32c) 



(33) 



and the corresponding averages for and are given 

by, 



coth 



(34) 



2A/2^2 



M7 



(35) 



From Eq. [26l putting Uq ^ we find that we can re- 
cover Eq. 1321 because in this limit a — > 1 and a' — > ag. 
However, it is puzzling to note that Eq. [35] cannot be 
recovered from Eq. I31bl by taking the limit of ojq 0. 
In order to recover the well known diffusive behaviour 
of the free particle, we must consider the system as it 
approaches equilibrium, viz. Eq. 1301 instead of the equi- 
librium values given by Eq. EH In Eq. [3D1 as 0, 
one of the exponents (7 — Q) which sets the relaxation 
time vanishes. Thus, we expand this term as a power 
series, i.e.,exp(-2(7- f])i) - l + {uj^/2j)t + 0{uj^). This 
term comes with the integration cosntant C2, which can 
be calculated as. 



_ (Dppfl + [2m^-/^ - 2m^-fujl + 2711^-^"^^ - m^f^wg]!?,, - 2raDpq[Vi^ + 7!^]) 

{{p^)t=o[^^ - 7^?] - {x^)t=om'^ujlP^ + 7^] - {qp + pq)t=omujl^) ,„„^ 



Inserting the power series expansion and the value of C2 
from Eq. [351 in Eq. [3D1 we recover Eqs. [Mj and [3S| ex- 
actly, in the limit of wq 0. Thus, it is clear that the 
equilibrium properties of a particle in a harmonic oscil- 
lator has an extra degree of freedom and hence the free 
particle results cannot be recovered in a straightforward 
manner. 



V. CONCLUSIONS 

We have considered a microscopic model of a parti- 
cle in a harmonic well, connected to a bath of oscillators 
and derived a master equation for the reduced density 
matrix of the system. We have calculated the diffusion 
constants in both the overdamped and the underdamped 



regimes and examined the positivity conditions for the 
reduced density matrix. The postivity condition is not 
satisfied below a breakdown temperature Tc, which de- 
pends on the dimcnsionless parameter wq/t- This is not 
very unexpected, since we neglected the transients for 
t 3> h/ksT. In the high temperature limit, we find that 
the diffusion constants become independent of the oscil- 
lator frequency (upto O (1/T)). The kinetic and poten- 
tial energy obtained from the steady state solutions are 
in good agreement with the thermodynamic results. We 
have also recovered the free particle results, in the limit 
of vanishing oscillator frequency. To recover the correct 
diffusive behaviour of the free particle, it is important 
to consider the approach to equilibrium, rather than the 
asymptotic equilibrium values. 
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